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Abstract 

Accurate spectroscopic constants and electrical properties of small molecules are determined by 
means of W4 and post-W4 theories. For a set of 28 first- and second-row diatomic molecules 
for which very accurate experimental spectroscopic constants are available, W4 theory affords 
near-spectroscopic or better predictions. Specifically, the root-mean-square deviations (RMSD) 
from experiment are 0.04 pm for the equilibrium bond distances (r e ), 1.03 cm -1 for the harmonic 
frequencies (u) e ), 0.20 cm -1 for the first anharmonicity constants (co e x e ), 0.10 cm -1 for the second 
anharmonicity constants (u e y e ), and 0.001 cm -1 for the vibration-rotation coupling constants (a e ). 
These RMSD imply 95% confidence intervals of about 0.1 pm for r e , 2.0 cm -1 for co e , 0.4 cm -1 for 
u} e x e , and 0.2 cm -1 for w e y e . We find that post-CCSD(T) contributions are essential to achieve 
such narrow confidence intervals for r e and u e , but have little effect on u; e x e and a e , and virtually 
none on u e y e - Higher-order connected triples, T3— (T), improve agreement with experiment for 
the hydride systems, but their inclusion (in the absence of T4) tends to worsen agreement with 
experiment for the nonhydride systems. Connected quadruple excitations, J4, have significant and 
systematic effects on r e , oj e , and u) e x e , in particular they universally increase r e (by up to 0.5 
pm), universally reduce ui e (by up to 32 cm -1 ), and universally increase co e x e (by up to 1 cm -1 ). 
Connected quintuple excitations, T5, are spectroscopically significant for oj e of the nonhydride 
systems, affecting oj e by up to 4 cm -1 . Diagonal Born-Oppenheimer corrections have systematic 
and spectroscopically-significant effects on r e and uj e of the hydride systems, universally increasing 
r e by 0.01-0.06 pm and decreasing uj e by 0.3-2.1 cm -1 . Obtaining r e and u e of the pathologically 
multireference BN and BeO systems with near-spectroscopic accuracy requires large basis sets in 
the core- valence CCSD(T) step and augmented basis sets in the valence post-CCSD(T) steps in W4 
theory. The triatomic molecules H2O, CO2, and O3 are also considered. The equilibrium geometries 
and harmonic frequencies (with the exception of the asymmetric stretch of O3) are obtained with 
near-spectroscopic accuracy at the W4 level. The asymmetric stretch of ozone represents a severe 
challenge to W4 theory, in particular the connected quadruple contribution converges very slowly 
with the basis set size. Finally, the importance of post-CCSD(T) correlation effects for electrical 
properties, namely dipole moments (/x), polarizabilities (a), and first hyperpolarizabilities (f3) is 
evaluated. 
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I. INTRODUCTION 



Since the early 1990s, there has been growing interest in calculating accurate ab initio 
anharmonic force fields for small molecules. There are two aspects to be considered: (i) 
the nuclear vibrational analysis, and (ii) the level of the electronic structure calculations. 
For diatomic molecules (i) is relatively straightforward using a Dunham analysis, and for 
rigid triatomic molecules, such as those considered in the present work, (i) entails no serious 
problems. As for (ii), for systems with mild to moderate nondynamical correlation effects, 
the accuracy of CCSD(T) basis set limit results is on the order of 5-10 cm -1 for harmonic 
frequencies^' and 0.1-0.5 pm for bond lengths. 1 ™ 3 To surpass this level of accuracy 
(or when considering molecules with more pronounced multireference character such as 
C2, BN, BeO, and O3) one should consider multireference methods®^ or post-CCSD(T) 
contributions.^^ Recently, Mintz et a/P have applied the multireference ccCA method^ 
to compute the potential energy curves of N 2 and C2; they obtained errors of ~0.06 pm for 
the bond lengths and ~2 cm^ 1 for the harmonic frequencies. 

The determination of very precise molecular energies (atomization energies or enthalpies 
of formation) has been one of the primary goals of composite ab initio methods. It has been 
shown that highly accurate thermochemical protocols such as HEAT^HEl anc [ w4 17 | 1S | are 
capable of sub-kJ/mol accuracy on average (see Ref P for a recent review). In particular, 
W4 theory obtains a root mean square deviation (RMSD) of 0.08 kcal/mol against a test 
set of 25 first- and second-row small molecules for which very accurate ATcT 20 atomization 
energies are available. This implies a 95% confidence interval of ~0.16 kcal/mol; the mean 
signed deviation (MSD) of just -0.01 kcal/mol suggests that W4 is free of systematic bias. 
Worst-case errors for problematic molecules are < 1 kJ/mol: for example, for F202,^C2p^ 
F 2 0, 21 and Og, 113 W4 is 0.09, 0.17, 0.18, and 0.23 kcal/mol away from the ATcT atomization 
energies (the ATcT value for C2 is available in Ref. and for the remaining species in RefP^j). 
In our continued work in this field we have shown that contributions from successively higher 
cluster expansion terms converge increasingly faster with the basis set, as they increasingly 
reflect nondynamical rather than dynamical correlation. 1 17 1 18 1 Indeed, the fact that W4 theory 
and related methods can be carried out at all at a realistic cost hinges on this behavior. 

In view of the success of W4 theory for thermochemical propertie d 17 * 18 * 21 ^ it is of interest 
to establish the reliability of W4 (and related) methods for spectroscopic properties based 
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on energy derivatives with respect to the nuclear coordinates (e.g. r e , u e , and u e x e ) and for 
electrical properties based on energy derivatives with respect to an external static electric 
field (e.g. ji, a, and f3). 

The present study considers a chemically diverse data set of 31 first- and second-row 
diatomic molecules as well as the triatomic molecules H2O, CO2, and O3. The chosen set, 
which includes radicals, polar systems, hydrides, and nonhydrides with single and multiple 
bonds, evidently spans a wide gamut from systems dominated by dynamical correlation (e.g. 
BH and H 2 0) to systems with pathological nondynamical correlation (e.g. C2, BN, BeO, 
and 3 ) and all shades in between. 

II. COMPUTATIONAL METHODS 
A. Electronic structure 

All the SCF, CCSD, and CCSD(T) calculations were carried out using MOLPRO 
2009.1^1 running on the Martin group Linux cluster at the Weizmann Institute. The post- 
CCSD(T) calculations were carried out using MRCC interfaced to the MOLPRO program 
suite. ^ The diagonal Born-Oppenheimer correction (DBOC) calculations were carried out 
using the CFOUR program system.^ All basis sets employed belong to the correlation 
consistent family of Dunning and co-workers The notation aug'-cc-pV(n-|-d)Z indicates 
the combination of regular cc-pVnZ on hydrogen and aug-cc-pV(n-|-d)Z on other elements. 

The computational protocols of W4, W4.2, and W4.3 theories have been specified and 
rationalized in great detail elsewhere. 1 17 1 18 1 In the present work the steps involved in W4 
theory are divided in the following manner: 

• W4[up to CCSD(T)j cnr ] represents the clamped-nuclei, non-relativistic CCSD(T) 
infinite basis set limit energy in the frozen-core approximation (in which the Is orbitals 
for first-row atoms and the Is, 2s, and 2p orbitals for second-row atoms are constrained 
to be doubly occupied in all configurations). The fc and nr subscripts stand for "frozen- 
core" and "non-relativistic". The following extrapolations are used for the HF, CCSD, 
and (T) contributions: 

— The ROHF-SCF contribution is extrapolated from the aug'-cc-pV(5+d)Z and 
aug'-cc-pV(6+d)Z basis sets using the Karton-Martin modification^ of Jensen's 
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extrapolation formula.' 33 ' 

— The RCCSD valence correlation energy is calculated from these same basis sets. 
Following the suggestion of Klopper ™ it is partitioned into singlet-coupled pair 
energies, triplet-coupled pair energies, and T± terms. The singlet-coupled and 
triplet-coupled pair energies are extrapolated using the A + B/L a two-point 
extrapolation formula (where L is the highest angular momentum present in the 
basis set) with 0:5=3 and 0-^=5, respectively, and the T\ term (which exhibits 
very weak basis set dependence) is set equal to that in the largest basis set. 

— The (T) valence correlation energy is extrapolated from the aug'-cc-pV(Q+d)Z 
and aug'-cc-pV(5+d)Z basis sets using the A + B/L 3 two-point extrapolation 
formula. For open-shell systems the Werner-Knowles-HampeP^ (a.k.a. 
MOLPRO) definition of the restricted open-shell CCSD(T) energy is employed 
throughout, rather than the original Watts-Gauss-BartlettF^ (a.k.a. ACES 
II/CFour) definition. 

• W4[up to CCSD(T)] in addition includes inner-shell correlation and scalar relativistic 
contributions. The former is extrapolated from RCCSD (T)/aug'-cc-pwCVTZ and 
RCCSD(T)/aug'-cc-pwCVQZ energies using the A + B/L 3 two-point extrapolation 
formula, and the latter (in the second-order Douglas-Kroll-Hess approximation^^) is 
obtained from the difference between non-relativistic CCSD(T)/aug'-pV(Q+d)Z and 
relativistic CCSD(T)/aug'-cc-pV(Q+d)Z-DK calculations. 1311 

• W4[up to CCSDT] additionally includes higher-order connected triples, T 3 — (T), 
valence correlation contribution extrapolated from the cc-pVDZ and cc-pVTZ basis 
sets using the A + B/L 3 two-point extrapolation formula. 

• In W4[up to CCSDTQ] the connected quadruples, T4, term is also included. The 
valence correlation (Q) and T 4 — (Q) contributions are calculated with the cc-pVTZ 
and cc-pVDZ basis sets, respectively. In Refs.^^ we found that scaling their sum by 
1.10 offers a very reliable (as well as fairly cost-effective) estimate of the basis set limit 
T 4 contribution. 

• Adding the connected quintuple, T 5 , valence correlation contribution calculated with 
the sp part of the cc-pVDZ basis set — denoted cc-pVDZ(no d) — results in full W4 
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theory. The T5 contribution converges very rapidly with the basis set as it primarily 
represents static correlation] 17 * 18 ! 

W4 represents an approximation to the relativistic basis-set limit CCSDTQ5 energy. The 
DBOC contributions, calculated at the ROHF/aug'-cc-pVTZ level of theory, are reported 
separately in Tables [I}fV] and are included in the final error statistics presented in these 
tables. 

For the smaller systems we also consider the post-W4 methods (W4.2 and W4.3) as 
defined in RefP^The changes in W4.2 and W4.3 relative to W4 are summarized as follows: 

• W4.2 theory in addition takes account of the T 3 — (T) correction to the core- valence 
contribution obtained using the cc-pwCVTZ basis set. 

• In W4.3 all the valence post-CCSD(T) corrections are additionally upgraded: the 
T 3 — (T) and (Q) corrections are extrapolated from the cc-pVTZ and cc-pVQZ basis 
sets, the T 4 — (Q) and T 5 corrections are calculated with the cc-pVTZ and cc-pVDZ 
basis sets, respectively, and the T 6 correction is calculated with the cc-pVDZ(no d) 
basis set. 



Finally, for the pathologically multireference systems in Table VI (C2, BN, and BeO) we 
also consider two additional extensions of W4 theory: 

• Replacing the regular correlation-consistent basis sets on electronegative atoms (N and 
O) with their augmented versions in all the valence post-CCSD(T) steps in W4 theory. 

• Extrapolating the CCSD(T) inner-shell correlation contribution from larger basis sets, 
namely the aug-cc-pwCVQZ and aug-cc-pwCV5Z basis sets. 



B. Vibrational analysis 

For each of the diatomic molecules, a 21-point potential energy curve was calculated using 
the above-mentioned W4 (and related) theories. The single point energy calculations are 
carried out at bond distances r x = rl xpt + x (x = 0, ±1, ±2, . . . , ±10 pm), where r expt " is the 
experimental equilibrium bond length. The spectroscopic constants (r e , u e , u e x e , u e y e , and 
a e ) are obtained from a Q th -g th degree Dunham fit of the potential energy curves.^ 
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For the triatomics CO2, H2O, and O3 grids of all points required to generate the 
nonvanishing quadratic, cubic, and quartic force constants in symmetry-adapted internal 
coordinates are calculated using the above-mentioned W4 (and related) theories. 

The convergence criteria throughout are tightened such that the SCF and coupled-cluster 
energies are converged to at least 1CT 12 a.u. 



C. Response properties 

The dipole moments (//), static polarizabilities (a), and static hyperpolarizabilities (/3) 
are calculated at the experimental structures via numerical differentiation of energies with 
respect to a uniform electric field. A field strength of /i=0.0025 a.u. was chosen as the 
basic grid size: the values reported were obtained by means of Richardson extrapolation^ 
from successive multiples of that grid size. The first-order (/i 2 ) and third-order (/3 ZZZ ) partial 
derivatives are obtained by means of a four point formula using static field strengths of ±/t 
and ±2h. The second-order derivatives (a xx , a yy , and a zz ) are obtained through a five point 
formula using field strengths of 0, ±h, and ±2h. The mixed third-order derivatives ((3jj Z , j = 
x,y) are computed using a six point formula with field strengths of ±2hz, ±2/i(z + j), and 
±2/i(2 -j). 



III. RESULTS AND DISCUSSION 



A. Multireference considerations 



In reference^ we proposed a number of energy-based diagnostics for the importance of 
nondynamical correlation that are specifically designed for thermochemical purposes. In 
particular, the %TAE e [(T)] diagnostic, i.e. the percentage of the total atomization energy 
accounted for by parenthetical connected triple excitations, proved to be a very efficient and 
cost-effective a priori diagnostic: 

Wmlnr/m , TAE e [CCSD(T)l — TAEJCCSD] 

%TAE e [(T)] = 100 x TAEe[CCSD(T)] ' « 

where TAE e [CCSD] and TAE e [CCSD(T)] represent, respectively, the non-relativistic, 
clamped-nuclei CCSD and CCSD(T) total atomization energies at the bottom of the well. 
Values of %TAE e [(T)] are interpreted as follows^ below 2% indicates systems dominated by 



dynamical correlation; 2-5% mild nondynamical correlation; 5-10% moderate nondynamical 
correlation; and in excess of 10% severe nondynamical correlation. 

Table S-I of the Supporting Information 46 provides the percentage of the non-relativistic, 
clamped-nuclei total atomization energy at the bottom of the well accounted for by SCF, 
(T) triples, post-CCSD(T), and T4 + T5 contributions. Table S-I also lists the coupled cluster 
T\ diagnostic,^ £> a diagnostic,^ and the largest CCSD T 2 amplitudes. However, the use of 
these diagnostics as a measure of multireference effects has been criticized by u d 17 * 21 ! and by 
others.HSlHSIln the context of the present work, for example, a 7i diagnostic of 0.03 is obtained 
for CS and SiO as well as for O3, although the latter is a notoriously multireference system. 
The V>\ diagnostic of O2 (0.01) is the same as, or even smaller than, that of the hydride 
systems, despite the former having almost an order of magnitude greater %TAE e [T4 + T5]. 
The largest T 2 cluster amplitudes are deceptively high for BH (0.11) and deceptively low for 
BeO (0.04). In practice, the %TAE[(T)] indicator proves to be a very useful predictor for 
the importance of post-CCSD(T) correlation effects (see relevant discussions in RefsJ^^j). 

As expected, the hydride systems considered in the present study exhibit very mild 
nondynamical correlation effects and can be regarded as dominated by dynamical correlation: 
60-80% of the atomization energy is accounted for at the SCF level, and only 0-2% by the 
(T) triples. The nonhydrides BF, BC1, A1F, A1C1, CO, SiO, SiS, and N2 are characterized by 
mild nondynamical correlation: at least 50% of the TAE e is accounted for at the SCF level, 
and 2<%TAE e [(T)]<5. The homonuclear diatomics (0 2 , Si 2 , P 2 , S 2 , and Cl 2 ), as well as 
CS, PN, SO, C1F, and BeO are characterized by moderate nondynamical correlation effects, 
20<%TAE e [SCF]<60 and 6<%TAE e [(T)]<8. Finally, C 2 , BN, and 3 are pathologically 
multireference: at the SCF level BN and O3 are unbound, and only 10% of the atomization 
energy of C 2 is accounted for, while the (T) triples account for as much as 13-19% of the 
binding energies of these systems. 

B. Equilibrium geometries (r e ) of diatomic molecules 

Table |l| compares the equilibrium bond distances obtained at the W4 (and related) 
levels with experimental bond lengths for our set of 10 monohydride and 18 nonhydride 
diatomic molecules. The equilibrium bond lengths determined at the clamped-nuclei, 
frozen-core, non-relativistic CCSD(T) basis set limit (W4[up to CCSD(T)/ Cjnr ] method) 
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almost universally overestimate the experimental bond distances (by ~0.3 pm on average). 
Inclusion of the core-valence and scalar relativistic corrections consistently reduces r e by 
amounts ranging from ~0.1 pm (for HF, OH, F 2 , and NH) to ~1.0 pm (for A1C1 and 
Si2). Consequently, at the clamped- nuclei, all-electron, relativistic CCSD(T) basis set limit 
(W4[up to CCSD(T)] method) r e are universally underestimated (by ~0.2 pm on average). 
The largest underestimations are seen for the nonhydrides, e.g. by 0.26, 0.28, 0.31, 0.32, 
0.33, 0.35, 0.38, 0.39, and 0.52 pm, for SiS, PN, C1F, Cl 2 , 2 , S 2 , F 2 , P 2 , and Si 2 , respectively. 
For the subset of 18 nonhydride systems the W4[up to CCSD(T)] method attains an RMSD 
of 0.27 pm. As expected, the performance for the subset of monohydride systems, which 
are mostly dominated by dynamical correlation (0.1<%TAE[(T)]<1.6), is much better with 
an RMSD of only 0.06 pm. For the entire set of 28 diatomics W4[up to CCSD(T)] theory 
attains an RMSD of 0.22 pm. 

Higher-order triple excitations, T 3 — (T), have a relatively small effect on the bond lengths: 
\r e \ changes by <0.1 pm (with the exception of P 2 and PN for which r e decreases by 0.2 pm). 
For the hydride systems the T 3 — (T) contribution consistently increases the bond lengths by 
0.01-0.06 pm thereby improving agreement with experiment. For the nonhydride systems 
the T 3 — (T) excitations have a mixed effect, bettering agreement with experiment for some 
while worsening it for others. For the entire set of 28 diatomics, the effect of including the 
T 3 — (T) contribution is manifested by a slight increase in the RMSD from 0.22 to 0.25 pm. 
Note that, with very few exceptions (namely BH, CH, and A1H), the W4[up to CCSDT] 
method still universally underestimates the experimental bond lengths. 

Connected quadruple excitations, T 4 , universally increase the bond distances, thereby 
substantially improving the agreement with experiment, particularly for the nonhydride 
systems. Upon inclusion of the T4 excitations the RMSD over the nonhydrides is reduced 
from 0.30 to 0.06 pm (!), where in most cases the bond lengths increase by > 0.2 pm (e.g. by 
0.28, 0.29, 0.30, 0.33, 0.37, 0.42, and 0.46 for SO, S 2 , Si 2 , 2 , PN, F 2 , and P 2 , respectively). 
For the monohydrides, the effect of the T 4 excitations is rather modest, the bond lengths 
increasing by only 0.00-0.03 pm and the RMSD being reduced from 0.04 to 0.03 pm. 

As expected, connected quintuple excitations have essentially no effect on the internuclear 
distances of the hydride systems. For the nonhydride systems the bond lengths generally 
increase by 0.00-0.05 pm upon inclusion of the T 5 excitations. 

At the full W4 level the RMSD from experiment are {0.03, 0.05, 0.04} pm for the hydrides, 
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nonhydrides, and the entire set, respectively. The largest deviations (underestimations) are 
0.07 pm (for P 2 , SiH, SiS, and C1F) and 0.09 pm (for Cl 2 ). 

The trends discussed above can be qualitatively understood by looking at how the various 
energy contributions (£j core - valence ) E T3 ~^\ E Ti , and E Ts ) change as a function of r. Tables 
SII-SV of the Supporting Information® give these energy contributions as a function of r. 
The said energy contributions (with the exception of _E T3 ~( T ) for BF and A1F; vide infra) 
show a nearly linear dependence on r in the interval \r — r e \ < 10 pm. Addition of a linear 
term (E « ar) to a simple harmonic potential results in an increase in r e if a < and a 
decrease in r e if a > 0. For all the systems, the connected quadruples energy contribution, 
E Ti , decreases with r (in a nearly linear manner: R 2 > 0.99) and thereby inclusion of the T4 
excitations universally lengthens the bond distances. The core-valence energy contribution, 
^core- valence ^ systematically increases with r (in a fairly linear fashion, R 2 > 0.98), thus 
universally shortening the bond lengths. For the hydride systems E T3 ~( T ' decreases linearly 
with r (R 2 > 0.98), therefore lengthening the equilibrium bond lengths. For the nonhydrides, 
the E Ta "( T ) energy contribution does not behave in a consistent manner: for some (namely 
Si2, CI2, BC1, C1F, A1C1, and SiS) it decreases with r therefore increasing r e , while for others 
it increases with r, thus decreasing r e . Interestingly, for the highly polar BF and A1F systems 
^Ta-(T) var i es quadratically with r in the scanned interval (R 2 > 0.998, see Table S-II of 
the Supporting Information®) and r e is not affected by the T 3 — (T) correction. For most of 
the systems E Ts decreases with r (exceptions are: BF, BC1, A1F, A1C1, CS, SiS, and SiO), 
therefore inclusion of T 5 excitations increases r e for all but the said exceptions. In the next 
section we shall see that the said energy contributions that increase r e decrease u e and vice 
versa — as expected since uj oc -. 

For a few systems we were able to obtain potential curves at the W4.2 and W4.3 levels 
(see Table [i]). Adding the T 3 — (T) correction to the core- valence contribution in W4.2 theory 
has little effect on the bond distances, namely r e is increased by 0.02 and 0.03 pm for F 2 
and P2, respectively. Likewise, computing all the valence post-CCSD(T) corrections with 
larger basis sets and adding a T$ correction in W4.3 theory increases the bond length of F 2 
by merely 0.03 pm. 

Finally, DBOC corrections universally increase the bond distances and, as expected, only 
have a perceptible effect on the hydride systems, something more pronounced for the lighter 
hydrides within each row. For example, the DBOC correction increases r e by 0.06, 0.04, 
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0.03, 0.02, and 0.01 pm for the first-row hydrides BH, CH, NH, OH, and HF, respectively, 
and by 0.05, 0.04, 0.03, 0.02, and 0.01 pm for the second-row hydrides A1H, SiH, PH, SH, 
and HC1, respectively. 



C. Harmonic frequencies (w e ) of diatomic molecules 

The harmonic frequencies at the W4 (and related) levels are given in Table [XXj together 
with the experimental values. At the clamped-nuclei, all-electron, relativistic CCSD(T) 
basis set limit (W4[up to CCSD(T)] method) the harmonic frequencies are universally 
overestimated by amounts ranging from ~1 cm -1 (for A1F and A1C1) up to 29 cm -1 (for 
O2). Not surprisingly, the smallest errors (of 1-8 cm -1 ) are observed for the systems for 



which %TAE e [(T)]<4 (see section jlTTA] ) . The RMSD of the W4[up to CCSD(T)] method 
are {5.1, 11.2, 9.8} cm -1 for the hydrides, nonhydrides, and the entire set, respectively. 

For most systems W4[up to CCSD(T)/ Cinr ] and W4[up to CCSD(T)] err on the same 
side of experiment, the exceptions are A1F, CO, BC1, CH, BF, and BH for which W4[up 
to CCSD(T)/ Cjnr ] underestimates experiment by 2.5-8.6 cm -1 . We note that with few 
exceptions (namely BH, A1H, BF, A1F, CH, and SiH) W4[up to CCSD(T) /Cjnr ] is closer 
to experiment than W4[up to CCSD(T)], as manifested in an RMSD for the entire set 
smaller by 2 cm -1 . 

Higher-order triple excitations, T3— (T), consistently reduce the harmonic frequencies of 
the monohydride compounds by 1.4-5.0 cm -1 , thus cutting the RMSD for the monohydrides 
by almost 50% from 5.1 cm" 1 (W4[up to CCSD(T)]) to 2.8 cm" 1 (W4[up to CCSDT]). For 
the nonhydrides, the higher-order triple excitations have a mixed effect, reducing u e for 
some systems (e.g., by ~1 cm' 1 for S12, BC1, and CI2) and increasing u e for others (most 
notably by ~10 cm -1 for N 2 and PN). In effect, the RMSD for the nonhydrides increases 
upon inclusion of the T3— (T) correction (from 11.2 cm -1 for W4[up to CCSD(T)] to 15.2 
cm" 1 for W4[up to CCSDT]). 

Connected quadruple excitations, T 4 , universally reduce u e by amounts ranging from 
0.0-5.0 cm -1 for the monohydride systems and 0.3-31.5 cm -1 for the nonhydride systems. 
The largest decreases, of over 10 cm -1 , are seen for P 2 , SO, F 2 , N 2 , PN, and O2 (specifically 
by 13.4, 16.1, 16.7, 20.4, 21.8, and 31.5 cm -1 , respectively). Inclusion of the T4 excitations 
reduces the RMSD over the set of hydrides from 2.8 cm -1 (W4[up to CCSDT]) to 0.6 
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cm 1 (W4[up to CCSDTQ]), and over the set of nonhydrides from 15.2 cm 1 (W4[up 
to CCSDT]) to 2.4 cm" 1 (W4[up to CCSDTQ]). Over the entire set of 28 diatomics the 
RMSD is reduced from 12.7 cm" 1 (W4[up to CCSDT]) to 2.0 cm" 1 (W4[up to CCSDTQ]), 
an improvement by nearly one order of magnitude (!). This clearly demonstrates that 
connected quadruple excitations are essential for obtaining the harmonic frequencies with 
near-spectroscopic accuracy. 

As expected, connected quintuple excitations, T 5 , have little or no effect for the hydrides: 
the biggest changes are seen for OH and HF for which u e is reduced by merely 0.1 cm -1 . 
For the nonhydrides T 5 excitations change u e by up to 4 cm -1 . The biggest changes are 
observed for N 2 , PN, and O2, namely, u e is reduced by 3.3, 3.5, and 4.0 cm -1 , respectively. 

Upon inclusion of the T5 excitations the RMSDs for the nonhydrides (and for the entire 
set) are reduced by 50%, thus at the W4 level RMSDs of {0.6, 1.2, 1.0} cm -1 are obtained for 
the hydrides, nonhydrides, and the entire set, respectively. The largest deviation being 2.1 
cm -1 for F 2 and P 2 . We note that these errors drop to 0.9 and 1.3 cm -1 , respectively, upon 
adding the T 3 — (T) correction to the core- valence contribution in W4.2 theory. The error 
for F 2 is further reduced to 0.5 cm -1 by upgrading the valence post-CCSD(T) corrections 
and adding a Tq correction in W4.3 theory. 

Finally, a word is due on the effect of the DBOC correction on the harmonic frequencies. 
For the monohydride systems the DBOC reduces u e across the board, again more appreciably 
so for the lighter hydrides within each row. Specifically, the DBOC reduces u e by 2.1, 1.8, 
1.4, 0.9, and 0.3 cm -1 for the first-row hydrides BH, CH, NH, OH, and HF, respectively, and 
by 1.1, 1.0, 0.9, 0.7, and 0.5 cm" 1 for the second-row hydrides A1H, SiH, PH, SH, and HC1, 
respectively. For the nonhydrides the DBOC negligibly increases the harmonic frequencies 
by <0.01 cm -1 , the biggest change of +0.08 cm -1 being found for 2 . 



D. First anharmonicity constant (u e x e ) of diatomic molecules 



Table HI lists the theoretical and experimental first anharmonicity constants. 
Monohydrides are well known to exhibit relatively strong anharmonicity effects: u e x e ranges 
between 30-90 cm -1 for the subset of monohydrides, and between 2-14 cm -1 for the subset 
of nonhydrides. 

The W4[up to CCSD(T)] method almost universally underestimates u e x e by amounts 
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ranging from 0.01 cm" 1 (BF) to 1.0 cm" 1 (NH). At this level RMSDs of {0.6, 0.3, 0.5} 
cm -1 are obtained for the hydrides, nonhydrides, and the entire set, respectively. We note 
that the W4[up to CCSD(T) f c ,nr] method shows very similar performance for most of the 
systems, the core-valence correction tends to slightly increase u e x e but overall offers little 
improvement. 

Inclusion of higher-order triple excitations improves agreement with experiment for the 
monohydride systems, and slightly worsens agreement with experiment for the nonhydride 
systems. For the monohydride systems T 3 — (T) excitations increase u e x e by up to 0.5 cm -1 
(for NH), whereas for the nonhydrides they tend to decrease u e x e by up to 0.3 cm -1 (for 
N 2 and F 2 ). The RMSD of the W4[up to CCSDT] method is 0.4 cm" 1 for both the hydride 
and nonhydride subsets. 

Inclusion of connected quadruple excitations in the W4[up to CCSDTQ] method 
universally increases u e x e by up to 1.0 cm -1 . Thus, improving agreement with experiment 
across the board, something more pronounced for the nonhydride systems. The largest 
improvements are seen for N 2 and 2 : the error for N 2 is reduced from 0.7 cm -1 to 0.2 
cm -1 , and that for 2 from 1.1 cm -1 to 0.1 cm -1 . 

Contributions from connected quintuple excitations are generally negligible (for most of 
the systems <0.01 cm" 1 ) the largest changes being +0.06, +0.11, and +0.15 cm" 1 for PN, 
N 2 , and 2 , respectively. 

For the nonhydride systems agreement between W4 and experiment is as good as one 
can hope for with an RMSD of 0.04 cm -1 . For the hydride systems the RMSD is almost 
one order of magnitude larger (0.32 cm" 1 ) but still well within the goal of "spectroscopic 
accuracy" . 

W4.2 and W4.3 theories offer no significant improvement over W4 theory. Likewise, 
DBOC corrections are negligible (the largest correction of -0.02 cm -1 is obtained for BH). 



E. Second anharmonicity constant (u} e y e ) of diatomic molecules 



Table IV reports the theoretical and experimental second anharmonicity constants. The 
second-order anharmonicity constants are typically 2-3 orders of magnitude smaller than 
the first-order anharmonicity constants. For the monohydrides ui e y e ranges from 0.1 cm -1 
(for SiH, PH, and SH) to 0.9 cm -1 (for HF). For the nonhydrides, which exhibit weaker 
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anharmonicity effects, u e y e are typically <0.01 cm 1 , and are given in Table S-VI of the 
Supporting Information. 

The W4[up to CCSD(T) /Cj „ r ] method attains near-zero RMSD of 0.003, 0.09, and 0.06 
cm -1 for the nonhydrides, hydrides, and the entire set, respectively. For both subsets 
core- valence, post-CCSD(T), and DBOC corrections have very little effect on the second 
anharmonicity constants. 

Figure [T] shows the normal distribution functions for the errors between the theoretical 
and experimental r e , co e , u e x e , and oj e y e . It illustrates that r e and u e are much more sensitive 
to post-CCSD(T) correlation effects, in particular to the T 4 excitations, than the anharmonic 
corrections. This was suggested long ago, based on CEPA^ calculations on the methane 
surface by Pulay and coworkers. 48 The Gaussians for the W4[up to CCSDT] and W4[up to 
CCSDTQ] methods are centered, respectively, around -0.17 and -0.04 pm for r e , around 8.1 
and 1.0 cm -1 for u e , around -0.3 and -0.1 cm -1 for u e x e , and around -0.03 and -0.04 cm -1 
for u e y e . 



F. Vibration-rotation coupling constant (a e ) of diatomic molecules 

Table [V] depicts the theoretical and experimental vibration-rotation coupling constants. 
The RMSD of the W4[up to CCSD(T) /C)W .], W4[up to CCSD(T)], W4[up to CCSDT], 
W4[up to CCSDTQ], and W4 methods are 0.0034, 0.0024, 0.0017, 0.0012, and 0.0012 cm" 1 , 
respectively. Hence it can be concluded that, while the a e are reasonably well converged at 
the clamped-nuclei, frozen-core, non-relativistic CCSD(T) basis set limit (W4[CCSD(T)/ Cjnr .] 
method), significant improvements can still be obtained beyond it. 



G. Difficult 8- valence electron systems: C2, BN, and BeO 

The low-lying 1 S + electronic states of the twelve-electron isoelectronic diatomics C2, 
BN, and BeO are known to have significant contributions from singly and doubly excited 
determinants, and thus represent an extremely challenging test for single reference electron 
correlation methods.^ Table VI lists the spectroscopic constants at the W4 (and related) 



levels together with the experimental values. The harmonic frequencies are severely 
overestimated at the W4[up to CCSD(T)] level, specifically by 17, 50, and 10 cm -1 for 
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C2, BN, and BeO, respectively. Higher-order triples, T 3 — (T), contributions increase the 
errors for C2 and BeO to about 20 cm -1 , but reduce the error for BN by 36 cm -1 (!) to 
14 cm -1 . (Interestingly, the same decrease of 36 cm -1 was computed nearly two decades 
ago by one of us using a double zeta quality basis set!®) As expected, connected quadruple 
excitations reduce u e by substantial amounts, namely 15.5, 6.6, and 12.6 cm -1 for C2, BN, 
and BeO, respectively. Note however, that for BN the T4 excitations have a lesser effect on 
u e than the T3— (T) excitations. T5 excitations further reduce the errors for C2 and BN by 
1.5 and 2.0 cm -1 , respectively, but increase the error for BeO by 7.2 cm -1 . Thus, at the W4 
level the harmonic frequencies of C2, BN, and BeO are overestimated by 2.6, 5.7, and 14.3 
cm -1 , respectively. 

We suspected that the large discrepancies seen for the ionic BN and BeO molecules are 
partly due to the difficulty of the unaugmented basis sets to describe the charge distribution 
in the valence post-CCSD(T) steps. Adding diffuse functions on the electronegative atoms 
(N and O) in all the valence post-CCSD(T) steps in W4 theory (see Computational Methods 
section and footnote b in Table |VI[ ) reduces the errors for BN and BeO to 2.1 and 6.1 
cm -1 , respectively. In addition, because of the fairly small core-valence gap in BeO, we 
also considered extrapolating the CCSD(T) inner-shell correlation contribution from larger 



basis sets (see Computational Methods section and footnote c in Table VI), which results 
in further reduction of the error for BeO to 1.2 cm -1 . (The reduction to 1.9 cm -1 for BN is 
not significant.) 



As mentioned in section III B the W4[up to CCSD(T)] method consistently 
underestimates the experimental bond lengths of the diatomic molecules in Table [T] (by 
~0.2 pm on average). Similarly, the bond distances of C2 and BeO are underestimated 
by 0.2 pm, but that of BN is underestimated by as much as 1.0 pm. Higher-order triples, 
T3— (T), have little effect on r e of C2 and BeO, but increase r e of BN by as much as 0.7 pm 
(!) narrowing the gap between theory and experiment to 0.3 pm. T4 contributions further 
reduce the discrepancy between theory and experiment: W4[up to CCSDTQ] underestimates 
the experimental bond lengths of C2, BN, and BeO by 0.05, 0.12, 0.11 pm, respectively. T 5 
excitations increase r e of C2 and BN by 0.03 and 0.05 pm, but decrease r e of BeO by 0.11 
pm. Thus, at the W4 level the C2 and BN bond lengths are in relatively good agreement 
with experiment (underestimating by 0.03 and 0.07 pm, respectively), but r e of BeO is 
underestimated by as much as 0.21 pm. Again, augmenting the basis sets on N and O in all 
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the valence post-CCSD(T) steps in W4 theory reduces the discrepancy for BeO to 0.08 pm, 
and for BN to 0.03 pm. Extrapolating the inner-shell correlation from RCCSD(T)/aug-cc- 
pwCVQZ and RCCSD(T)/aug-cc-pwCV5Z calculations results in further reduction of the 
error for BN and BeO to merely 0.02 pm. 

Turning our attention to the first anharmonicity constants, for C2 errors of 1.1 and 0.1 
cm -1 are obtained at the W4[up to CCSD(T)] and W4 levels, where the biggest improvement 
comes from the T 4 excitations (see Table VI). For BN errors of 5.5 and 1.2 cm -1 are obtained 



at the W4[up to CCSD(T)] and W4 levels, but here the improvement essentially stems from 
the T 3 — (T) contribution. For BeO errors of 0.3 and 0.5 cm -1 are obtained at the W4[up to 
CCSD(T)] and W4 levels. We note that by adding diffuse functions on N and O in all the 
post-CCSD(T) steps in W4 theory the errors for BN and BeO are further reduced to 0.85 
and 0.05 cm -1 , respectively. 



H. Triatomic systems 

The W4 equilibrium geometries, harmonic frequencies, and fundamental frequencies of 



H 2 0, CO2, and O3 are given in Table VII together with the available experimental values. 
The W4 equilibrium geometries are in close agreement with the experimental ones. The 
equilibrium bond distances of CO2 and O3 are both 0.02 pm lower than experiment, and the 
bond angle of O3 is 0.07° higher than that of Barbe et al. 49 For H2O the W4 bond distance 
(95.76 pm) is bracketed in between the experimental value of Benedict et al™ (95.72 pm) 
and the more recent value of Jensen^ (95.84 pm), and the W4 bond angle (104.50°) is in 
between the said experimental values (104.52°®! and 104.44°^, respectively). The Born- 
Oppenheimer correction is expected to be negligible for CO2 and O3, and is on the order of 
0.003 pm and 0.02° for the bond distance and bond angle of water, respectively.^ We note 
that our W4 r e and 8 e for water differ by 0.02 pm and 0.02°, respectively, from the ab initio 
Born-Oppenheimer geometry reported in RefP^l. 

The W4 harmonic frequencies for CO2 and H2O are in relatively good agreement 
with experiment. In particular, the symmetric, bending, and antisymmetric vibrational 
frequencies of CO2 underestimate the experimental values by 1.9, 2.4, and 0.1 cm -1 , 
respectively, and those of H2O overestimate experiment!^ by 2.8, 0.5, and 3.3 cm -1 , 
respectively. 
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Ozone represents an extremely challenging system for both single reference and 
multireference methods (see, e.g^^ and references therein). The theoretical symmetric 
stretch converges smoothly with the level of theory, i.e. it overestimates the experimental 
value of Barbe et alF^ by 20.3, 3.7, and 0.6 cm -1 , respectively, at the W4[up to CCSDT], 
W4[up to CCSDTQ], and W4 levels. The harmonic bending frequency converges more 
rapidly, in particular it overestimates the experimental one by 2.9, 0.5, and 0.2 cm" 1 
at the same levels, respectively. However, the asymmetric stretch (u^), which exhibits 
severe nondynamical correlation character^^, behaves in a less systematic manner. It 
underestimates the experimental value by 1.6, 17.1, and 20.2 cm -1 , respectively, at the 
said sequence of levels (where the good agreement between the W4[up to CCSDT] value 
and experiment is clearly fortuitous). It is likely that the counterintuitive deterioration 
in agreement with experiment upon inclusion of the connected quadruple excitations is 
due to basis set deficiencies. The hypersensitivity of u 3 to the treatment of connected 
quadruple excitations has been noted in RefP^ Determination of oo^ requires evaluation of 
the energy at only one point (in C s symmetry) aside from the equilibrium C% v structure. 
We were able to obtain the CCSDTQ/cc-pVTZ(no / Id) energies at these two points — cc- 
pVTZ(no f Id) denotes the sp part of the cc-pVTZ basis set combined with the d function 
from the cc-pVDZ basis set. Replacing the T4 term in W4 theory with T4/cc-pVTZ(no 
/ Id) results in an asymmetric stretch of 1078.1 cm -1 , i.e. the error from experiment is 
reduced by more than 50%. Additionally upgrading the basis set for the parenthetical 
connected quadruple excitations to cc-pVTZ {i.e., T 4 =CCSDT(Q)/cc-pVTZ-CCSDT/cc- 
pVTZ+CCSDTQ/cc-pVTZ(no / ld)-CCSDT(Q)/cc-pVTZ(no / Id)) results in a minor 
lowering of the asymmetric stretch by ~1 cm -1 , and a further improvement in the basis 
set to the cc-pVQZ(no g If) basis set has a similar effect. Unfortunately, carrying out 
the fully iterative CCSDTQ/cc-pVTZ calculations is beyond the capabilities of our current 
computational resources, these calculations involving 5xl0 9 and 10.5 xlO 9 amplitudes in Civ 
and C s symmetries, respectively. 

It is appropriate to note, at this point, that — in spite of the difficulties the potential 
energy surface of ozone presents for W4 theory - W4 performs very well for the 
thermochemistry of O3. In particular, the post-CCSD(T) components are practically 
converged at the W4 level, as evident from very small differences between W4 and W4.3 
theories. The post-CCSD(T) contributions to the total atomization energy of ozone in W4 
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and W4.3 theories are: f 3 -(T) = -1.34 and -1.37; T 4 = 3.81 and 3.84; and f 5 = 0.41 and 
0.43.^ In addition, W4.3 includes a Tq contribution of 0.05 kcal/mol. This suggests that 
while basis set convergence is faster for second-order properties (e.g. frequencies) than for 
the energy, n-particle treatment convergence is arguably slower. 



I. Electrical properties 

Fairly accurate experimental dipole moments are available for some of the molecules 
considered in the present work; however, for the higher-order properties (a and 0) accurate 
experimental data are much more scarce. Furthermore, direct comparison between theory 
and experiment is hampered by the fact that the theoretical values correspond to the static 
limit, whereas the experimental values are generally frequency-dependent. Thus, we will 
not attempt to quantify the reliability of W4 theory against experimental data; our intent 
is rather to investigate the effects of core-valence and post-CCSD(T) contributions on the 
electrical properties. 

Table VIII| gives the theoretical and experimental dipole moments. Core-valence and 



post-CCSD(T) corrections have no noticeable effect on systems dominated by dynamical 
correlation effects, and relatively little effects on the strongly multireference systems. The 
dipole moments at the frozen-core, non-relativistic CCSD(T) basis set limit (W4[up to 
CCSD(T)j c „ r ] method) agree very well with the full W4 values (RMSD=0.005 a.u., and 
the largest deviation of 0.01 a.u. is seen for BeO). At the W4[up to CCSD(T)] level, the 
RMSD from W4 is 0.002 a.u., while the largest error of 0.005 a.u. is found for O3. Adding 
the higher order triple excitations increases the errors for the nonhydrides: in particular, 
a dramatic increase is seen for BN (from -0.003 to 0.05 a.u.). We note that for BN, the 
T 3 — (T) excitations raise /1 by 0.05 a.u. and the T 4 excitations reduce it by 0.04 a.u. At the 
W4[up to CCSDTQ] level, the largest errors of -0.014 and +0.012 a.u. are seen for BeO and 
BN, respectively, and the RMSD from W4 is 0.06 a.u. 



Table IX shows the theoretical polarizabilities together with the available experimental 
data. The RMSDs from the full W4 values are 2.6, 2.9, 0.5, and 0.2 a.u. for the W4[up to 
CCSD(T) /c , nr ], W4[up to CCSD(T)], W4[up to CCSDT], and W4[up to CCSDTQ] methods, 
respectively. However, for most of the systems in Table [lX| the polarizabilities at the W4[up 
to CCSD(T) f c ,nr] level agree to within 0.1 a.u. with the W4 values. The main exceptions are 
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a zz of BeO, and a zz and a xx of C2, which are overestimated by as much as 1.5, 5.4, and 13.0 
a.u., respectively. Upon excluding C 2 and BeO the RMSD of the W4[up to CCSD(T)j cnr ] 
method drops to 0.13 a.u. The W4[up to CCSD(T) /C;nr ] and W4[up to CCSD(T)] methods 
show very similar performance. Inclusion of T3— (T) excitations dramatically reduces the 
said errors for C2 to 0.1 and 1.0 a.u., respectively, and the RMSD for the entire set drops to 
0.5 a.u. Nevertheless, at this level errors of -1.0 and 2.1 a.u. are still seen for a xx and a zz of 
BeO, respectively. These are reduced to -0.3 and 1.0 a.u., respectively, by inclusion of the 
T 4 excitations in the W4[up to CCSDTQ] method. 

Table [X] lists the theoretical W4 (and related) hyperpolarizabilities. Similarly to 
the polarizabilities, for most systems, the first hyperpolarizabilities at the W4[up to 
CCSD(T)/ Cjnr ] level agree with the W4 values to within 0.1 a.u. The main exceptions 
are O3 {(3 yy z) and CS ((3 ZZZ ); excluding these, the W4[up to CCSD(T)/ Cjnr ] method attains 
an RMSD of 0.3 a.u. from the W4 values. Inclusion of the core-valence and relativistic 
contributions (in the W4[up to CCSD(T)] method) reduces the said RMSD to 0.1 a.u., and 
inclusion of post-CCSD(T) excitations offers little improvement. Turning our attention to 
the more problematic systems, we note that at the W4[up to CCSD(T) / c , nr ] level f3 zzz of CS 
is ~3 a.u. removed from the W4 value and that core-valence and post-CCSD(T) corrections 
do not significantly affect this result. The W4[up to CCSD(T)/ Ci „ r ] method gives the wrong 
sign for the (3 yyz component of O3, i.e. it overestimates the W4 value by as much as ~9 
a.u. Inclusion of T 3 — (T) excitations reduces this error to ~2 a.u. and further inclusion of 
T 4 excitations reduces the error to -0.2 a.u. 

J. Conclusions 

For a chemically diverse set of 28 first- and second-row diatomic molecules for which 
very accurate experimental spectroscopic constants are available W4 theory is capable of 
spectroscopically accurate predictions. Specifically, the RMSD from experiment are 0.04 
pm for r e , 1.03 cm -1 for u e , 0.20 cm" 1 for u e x e , 0.10 cm -1 for u e y e , and 0.001 cm -1 for a e . 

Core-valence and post-CCSD(T) contributions have little effect on u e x e and a e , and 
especially on u e y e . As for r e and co e , higher-order triple excitations, T 3 — (T), improve 
agreement with experiment for the monohydride systems, but their inclusion (in the absence 
of T4) generally worsens agreement with experiment for the nonhydride diatomics. For 
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instance, upon inclusion of the T 3 — (T) corrections the RMSD in co e reduces from 5.1 to 2.8 
cm -1 for the subset of monohydride systems, and increases from 11.2 to 15.2 cm -1 for the 
subset of nonhydride systems. 

Connected quadruple excitations, T 4 , are essential to obtain r e and u e with RMSD of 
~0.05 pm and ~2 cm -1 , respectively. For example, inclusion of the T 4 excitations reduces 
the RMSD in oo e from 2.8 to 0.6 cm -1 for the subset of monohydrides, and from 15.2 to 2.4 
cm -1 for the subset of nonhydrides. The RMSD in r e are reduced from 0.04 to 0.03 pm 
for the hydrides subset, and from 0.30 to 0.06 pm for the nonhydrides subset. Moreover, 
inclusion of T 4 excitations results in systematic trends: namely r e universally increase (by 
up to 0.5 pm), u e universally decrease (by up to 31.5 cm -1 ), and u e x e universally increase 
(by up to 1.0 cm -1 ). 

Connected quintuple excitations, T 5 , are of minor importance for r e and u e x e , but are 
spectroscopically significant for oo e . Inclusion of a CCSDTQ5/cc-pVDZ(no d) correction 
term cuts the RMSD in the harmonic frequencies for our set of 28 diatomic molecules by 
50% from 2.0 to 1.0 cm -1 , practically all of the improvement being due to the nonhydride 
systems. 

DBOC corrections result in systematic trends and are of importance mainly for r e and 
uj e of the hydride systems, something more pronounced for the lighter hydrides within each 
row. For example, for our set of monohydrides DBOC universally increase r e by 0.06-0.01 
pm and universally decrease u e by 2.1-0.3 cm -1 . 

For the pathologically multireference systems (BeO and BN) obtaining near-spectroscopic 
accuracy requires: (i) augmenting the correlation-consistent basis sets on N and O in the 
valence post-CCSD(T) steps in W4 theory, and (ii) extrapolating the CCSD(T) inner-shell 
contribution in W4 theory from the aug-cc-pwCV{Q,5}Z basis set pair. 

The present study also considers water, carbon dioxide, and ozone. The experimental 
geometries and harmonic frequencies (with the exception of the asymmetric mode of ozone) 
are reproduced relatively accurately at the W4 level. However, the asymmetric stretch of 
ozone, which exhibits severe nondynamical character, is underestimated by as much as 20 
cm -1 at the W4 level. This deviation is reduced by more than 50% when the connected 
quadruples, T 4 , term in W4 theory is calculated with the cc-pVTZ(no / Id) basis set. 

We also calculated the dipole moments (//), polarizabilities (a), and first hyperpolariz- 
abilities (0). We find that: (i) excluding the pathologically multireference systems, the post- 
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CCSD(T) contributions are of little importance for these electrical properties; (ii) the core- 
valence and DBOC contributions are also of little importance in most cases; and (iii) that the 
sensitivity to post-CCSD(T) effects, particularly for systems with significant nondynamical 
correlation, increases with increasing order of derivative with respect to the static electric 
field, i.e., in the order: \i < a < (3. 
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TABLE I: Equilibrium bond distances (r e , pm) of diatomic molecules obtained from a Dunham 
analysis of potential energy curves computed at the W4 and related levels. 
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0.04 


0.03 


0.03 
















error statistics for nonhydridcs (18 sysl 


cms) 6 










MSD 


0.30 


-0.23 


-0.26 


-0.05 


-0.04 












MAD 


0.36 


0.23 


0.26 


0.05 


0.04 












RMSD 


0.43 


0.27 


0.30 


0.06 


0.05 
















error statistics for everything 


(28 systems) 6 










MSD 


0.27 


-0.16 


-0.17 


-0.04 


-0.03 












MAD 


0.30 


0.17 


0.18 


0.05 


0.04 












RMSD 


0.33 


0.22 


0.25 


0.05 


0.04 













a DBOC correction at the HF/AVTZ level. 'Mean signed deviations (MSD), mean absolute deviation (MAD), and root mean 
squared deviation (RMSD) from experiment, note that the DBOC corrections (cfr. footnote a) are included in the error 
statistics. 
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TABLE II: Harmonic frequencies (u e , cm 1 ) of diatomic molecules obtained from a Dunham 
analysis of potential energy curves computed at the W4 and related levels. 





W4[up to 


W4 


W4 


W4 


W4 


W4.2 


W4.3 


ADBOC 


Expt. 6 




CCSD(T) 


[up to 


[up to 


[up to 














C»r] 


CCSD(T)] CCSDT] CCSDTQ] 












BH 


2360.24 


23 1 0. jy 


Z369.13 


2309.08 


2369.08 


2369.10 


2368.69 


-2.08 


2366.7296 


CH 


2858.08 


ogee no 

2500. yy 


2563.0 1 


2502.03 


2862.51 






-1.83 


2860.7508 


NH 


3286.50 


3292.03 


O O O *7 OO 

328 1. 33 


3284.65 


3284.62 






-1.43 


3282.7200 


OH 


3742.91 


Q7/ie pro 

3 1 40.08 


3/44.1U 


3 f 3y.oo 


3739.48 






-0.90 


3737.7610 


HF 


4141.64 


4143.57 


4143.71 


/] 1 Q O 07 

4138.5 / 


4138.77 


4138.79 


4138.51 


-0.31 


4138.3850 


A1H 


1688.45 


1685.17 


1682.54 


1682.14 


1682.14 






-1.10 


1682.3747 


SiH 


2048.29 


OH/1 "7 (V7 

2U4f X)( 


2043.63 


zU4z.yi 


2042.90 






-1.04 


2042.5229 


PH 


2370.75 


OQ*7A "70 

23 (0. <y 


2360.80 


ZO04.41 


2364.41 






-0.91 


2363.7740 


SH 


2700.78 


OT'Al A1 

2/01.91 


2098.04 


n/OAf ^n 
2oy6.42 


2696.40 






-0.73 


2696.2475 


HC1 


2994.17 


2995.65 


OOO/I 9Q 

zyy4.zo 


nnoi £A 

2yyi.ou 


2991.48 


2991.60 


2990.94 


-0.53 


2990.9248 


N 2 


2362.34 


Z61 i.yi 


noon 07 


2301.53 


2358.54 


2357.89 


2358.21 


0.04 


2358.5700 


o 2 


1605.48 


1 end 00 
1609.22 


1617.24 


1585.76 


1581.78 






0.08 


1580.1610 


F 2 


931.17 


y32.34 


936.12 


919.44 


918.96 


917.82 


917.44 


0.02 


916.9290 


Si 2 


516.41 


517.67 


516.23 


511.23 


510.78 






0.00 


510.9800 


P 2 


787.79 


•7AA OO 

Z92.23 


f y ; .02 


784.26 


782.87 


782.10 




0.00 


780.7700 


s 2 


734.42 


735.99 


736.93 


728.16 


727.48 






0.00 


725.7102 


Cl 2 


564.79 


ODD. DO 


564.48 


561.28 


561.24 


560.94 




0.00 


559.7510 


BF 


1394.96 


14U0.U1 


1404.45 


1403.02 


1403.25 






0.04 


1402.1587 


CO 


2166.60 




2178.86 


2170.22 


2170.37 


2170.16 


2170.35 


0.05 


2169.7559 


BC1 


836.91 


843 27 


841.90 


840.74 


840.93 






0.03 


840.2947 


CS 


1285.98 


1292.50 


1293.26 


1284.95 


1285.24 


1285.45 




0.02 


1285.1546 


SiO 


1241.66 


1248.53 


1250.98 


1241.56 


1242.79 






0.01 


1241.5439 


PN 


1344.72 


1352.54 


1363.39 


1341.62 


1338.13 






0.02 


1336.9480 


so 


1160.68 


1163.56 


1168.78 


1152.64 


1151.10 






0.00 


1150.7913 


A1F 


799.79 


803.09 


802.84 


802.36 


802.44 






0.00 


802.3245 


C1F 


791.38 


792.11 


791.38 


784.96 


784.81 






0.00 


783.4534 


A1C1 


481.75 


482.80 


482.42 


482.11 


482.12 






0.00 


481.7747 


SiS 


752.61 


755.69 


755.59 


750.30 


751.09 






0.00 


749.6456 






error 


statistics for hydrides (10 systems) c 








MSD 


1.88 


4.66 


1.96 


-0.10 


-0.13 










MAD 


4.49 


4.66 


2.15 


0.48 


0.46 










RMSD 


4.81 


5.09 


2.80 


0.61 


0.60 










error statistics for nonhydrides (18 systems) c 


MSD 


4.61 


9.11 


11.57 


1.67 


0.97 










MAD 


6.41 


9.11 


11.57 


1.69 


1.00 










RMSD 


8.67 


11.22 


15.17 


2.35 


1.21 














error statistics for everything (28 syst 


ems) c 








MSD 


3.64 


7.52 


8.14 


1.04 


0.58 










MAD 


5.73 


7.52 


8.21 


1.26 


0.80 










RMSD 


7.79 


9.78 


12.69 


1.99 


1.03 











a DBOC correction at the HF/AVTZ level. b From RefEHl. c See footnote b of TablefE] 
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TABLE III: First-order anharmonic corrections {u e x e , cm 1 ) of diatomic molecules obtained from 
a Dunham analysis of potential energy curves computed at the W4 and related levels. 





W4[up to 


W4 


W4 


W4 


W4 


W4.2 


W4.3 


ADBOO 


Expt.k 






[up to 


[up to 


[up to 














fc,nr] 


CCSD(T)] CCSDT] CCSDTQ] 












BH 


48.92 


49.16 


49.33 


49.33 


49.33 


49.34 


49.33 


-0.02 


49.33983 


CH 


63.66 


63.89 


64.19 


64.26 


64.26 






-0.01 


64.4387 


NH 


77.93 


78.07 


78.53 


78.74 


78.74 






0.01 


79.04 


OH 


84.38 


84.48 


84.72 


85.01 


85.02 






0.00 


84.8813 


HF 


89.73 


89.83 


89.87 


90.09 


90.10 


90.10 


90.13 


0.01 


89.9432 


A1H 


28.75 


28.45 


28.59 


28.61 


28.61 






-0.01 


29.05098 


SiH 


35.67 


35.59 


35.84 


35.89 


35.89 






-0.01 


36.0552 


PH 


43.10 


43.09 


43.46 


43.58 


43.58 






0.00 


43.907 


SH 


47.98 


48.01 


48.25 


48.40 


48.40 






0.00 


48.742 


HC1 


51.89 


51.92 


52.02 


52.18 


52.18 


52.18 


52.16 


0.00 


52.8 


N 2 


13.86 


13.90 


13.64 


14.14 


14.24 


14.26 


14.24 


0.00 


14.324 


O2 


11.03 


11.07 


10.85 


11.81 


11.96 






0.00 


11.95127 


F 2 


11.45 


11.41 


11.07 


11.40 


11.40 


11.45 


11.44 


0.00 


11.3221 


S12 


1.93 


1.94 


1.92 


2.02 


2.03 






0.00 


2.02 


P2 


2.72 


2.71 


2.68 


2.83 


2.84 


2.86 




0.00 


2.835 


S2 


2.73 


2.73 


2.70 


2.85 


2.86 






0.00 


2.8582 


CI2 


2.61 


2.59 


2.62 


2.67 


2.67 


2.68 




0.00 


2.69427 


BF 


11.71 


11.81 


11.85 


11.86 


11.86 






0.00 


11.82106 


CO 


13.11 


13.17 


13.12 


13.33 


13.30 


13.30 


13.27 


0.00 


13.28803 


BC1 


5.37 


5.42 


5.44 


5.47 


5.47 






0.00 


5.4917 


cs 


6.32 


6.37 


6.36 


6.54 


6.52 


6.51 




0.00 


6.502605 


SiO 


5.84 


5.90 


5.83 


6.02 


5.99 






0.00 


5.97437 


PN 


6.56 


6.57 


6.44 


6.76 


6.81 






0.00 


6.8958 


SO 


6.16 


6.17 


6.09 


6.47 


6.51 






0.00 


6.4096 


A1F 


4.78 


4.86 


4.86 


4.86 


4.86 






0.00 


4.849915 


C1F 


4.86 


4.85 


4.84 


4.95 


4.95 






0.00 


4.9487 


A1C1 


2.08 


2.10 


2.10 


2.11 


2.11 






0.00 


2.101811 


SiS 


2.52 


2.53 


2.52 


2.61 


2.59 






0.00 


2.58623 


error statistics for hydrides (10 systems) c 


MSD 


-0.62 


-0.57 


-0.34 


-0.21 


-0.21 










MAD 


0.62 


0.57 


0.34 


0.27 


0.27 










RMSD 


0.68 


0.64 


0.41 


0.32 


0.32 












error statistics for nonhydrides (18 syst 


ams) c 








MSD 


-0.18 


-0.15 


-0.22 


-0.01 


0.01 










MAD 


0.19 


0.17 


0.22 


0.05 


0.03 










RMSD 


0.28 


0.26 


0.35 


0.07 


0.04 










error statistics for everything (28 systems) c 


MSD 


-0.34 


-0.30 


-0.26 


-0.08 


-0.07 










MAD 


0.35 


0.31 


0.27 


0.13 


0.12 










RMSD 


0.48 


0.45 


0.39 


0.21 


0.20 











a DBOC correction at the HF/AVTZ level. b From ReiW. c See footnote b of Table[f] 
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TABLE IV: Second-order anharmonic corrections (uj e y e , cm 1 ) of diatomic molecules obtained 
from a Dunham analysis of potential energy curves computed at the W4 and related levels. 





W4[up to 


W4 


W4 


W4 


W4 


W4.2 


W4.3 


ADBOC" 


Expt. 6 




CCSD(T) 


[up to 


[up to 


[up to 














fc,nr] 


CCSD(T)] CCSDT] CCSDTQ] 












BH 


0.384 


0.369 


0.359 


0.360 


0.360 


0.359 


0.367 


0.000 


0.362 


CH 


0.322 


0.314 


0.301 


0.298 


0.300 






-0.002 


0.3634 


NH 


0.290 


0.138 


0.113 


0.096 


0.096 






0.001 


0.367 


OH 


0.654 


0.518 


0.515 


0.507 


0.507 






0.003 


0.5409 


HF 


1.064 


1.052 


1.063 


1.056 


1.056 


1.056 


1.052 


0.002 


0.92449 


A1H 


0.238 


-0.030 


-0.031 


-0.032 


-0.032 






-0.001 


0.24762 


SiH 


0.149 


0.104 


0.095 


0.092 


0.092 






-0.002 


0.1254 


PH 


0.029 


0.012 


-0.007 


-0.015 


-0.014 






0.001 


0.1059 


SH 


0.065 


0.061 


0.050 


0.040 


0.040 






0.000 


0.1124 


HC1 


0.021 


-0.018 


-0.020 


-0.027 


-0.027 


-0.028 


-0.028 


0.001 


0.21803 


error statistics for hydrides (10 systems) 


MSD 


-0.015 


-0.084 


-0.093 


-0.099 


-0.099 










MAD 


0.075 


0.112 


0.121 


0.126 


0.125 










RMSD 


0.094 


0.147 


0.155 


0.160 


0.160 










error statistics for everything (28 systems) c 


MSD 


-0.005 


-0.032 


-0.034 


-0.038 


-0.038 










MAD 


0.030 


0.045 


0.049 


0.051 


0.051 










RMSD 


0.061 


0.095 


0.100 


0.103 


0.103 











a DBOC correction at the HF/AVTZ level. b From RefPO. c See footnote b of Table[f] 
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TABLE V: Vibration-rotation coupling constants (a e , cm 1 ) of diatomic molecules obtained from 
a Dunham analysis of potential energy curves computed at the W4 and related levels. 



MSD 
MAD 
RMSD 



W4[upto W4 W4 W4 

CCSD(T) [up to [up to [up to 
/Cjnr ] CCSD(T)] CCSDT] CCSDTQ] 



Wl 



W4.2 W4.3 ADBOC" 



Expt. 1 



BH 





41852 





42106 





42207 





42205 





42205 


0.42207 


0.42216 


-0.00027 


0.421565 


CH 





53115 





53345 





53503 





53543 





53544 






-0.00026 


0.53654 


NH 





64301 





64513 





64754 





64869 





64871 






-0.00018 


0.65038 


OH 





71515 





71697 





71830 





72005 





72009 






-0.00016 


0.7242 


HF 





78693 





78854 





78893 





79043 





79047 


0.79048 


0.79048 


-0.00002 


0.7933704 


A1H 





18372 





18561 





18625 





18636 





18636 






-0.00010 


0.1870527 


SiH 





21455 





21641 





21741 





21761 





21761 






-0.00006 


0.21814 


PH 





24869 





25042 





25180 





25223 





25223 






-0.00005 


0.25339 


SH 





27622 





27774 





27865 





27922 





27923 






-0.00005 


0.2799 


HC1 





30344 





30488 





30532 





30594 





30595 


0.30593 


0.30588 


-0.00005 


0.3069985 


N 2 





01697 





01703 





01682 





01722 





01729 


0.01731 


0.01729 


0.00000 


0.017318 


o 2 





01512 





01518 





01501 





01578 





01590 






0.00000 


0.0159305 


F 2 





01244 





01244 





01225 





01259 





01260 


0.01264 


0.01263 


0.00000 


0.0125952 


Si 2 





00128 





00130 





00130 





00134 





00134 






0.00000 


0.00135 


P 2 





00142 





00143 





00141 





00146 





00147 


0.00148 




0.00000 


0.00149 


S 2 





00153 





00155 





00154 





00159 





00159 






0.00000 


0.00159754 


Cl 2 





00145 





00145 





00145 





00147 





00147 


0.00148 




0.00000 


0.001516 


BF 





01888 





01901 





01905 





01906 





01906 






0.00000 


0.01904848 


CO 





01733 





01741 





01737 





01754 





01752 


0.01752 


0.01750 


0.00000 


0.01750513 


BC1 





00667 





00674 





00677 





00679 





00679 






0.00000 


0.0068124 


cs 





00580 





00585 





00585 





00594 





00593 


0.00593 




0.00000 


0.00591835 


SiO 





00494 





00499 





00496 





00505 





00504 






0.00000 


0.00503784 


PN 





00534 





00537 





00528 





00547 





00550 






0.00000 


0.0055337 


SO 





00557 





00560 





00556 





00575 





00577 






0.00000 


0.0057508 


A1F 





00490 





00499 





00500 





00500 





00500 






0.00000 


0.004984261 


C1F 





00425 





00426 





00426 





00433 





00433 






0.00000 


0.0043385 


A1C1 





00159 





00161 





00161 





00162 





00162 






0.00000 


0.001611082 


SiS 





00143 





00145 





00145 





00148 





00147 






0.00000 


0.00147313 



MSD 
MAD 
RMSD 



MSD 
MAD 
RMSD 



error statistics for hydrides (10 systems) c 



-0.00514 -0.00325 -0.00214 -0.00147 -0.00146 
0.00514 0.00325 0.00219 0.00152 0.00150 
0.00549 0.00379 0.00279 0.00192 0.00190 



error statistics for nonhydrides (18 systems) c 



-0.00016 -0.00012 -0.00016 -0.00002 -0.00001 
0.00016 0.00012 0.00016 0.00003 0.00002 
0.00024 0.00021 0.00028 0.00005 0.00002 



error statistics for everything (28 systems) c 



-0.00194 -0.00124 -0.00087 -0.00054 -0.00053 
0.00194 0.00124 0.00089 0.00056 0.00055 
0.00341 0.00236 0.00174 0.00119 0.00118 



a DBOC correction at the HF/AVTZ level. b From RefPD. c See footnote b of Table[f] 
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TABLE VI: Spectroscopic constants (r e in pm, all the rest in cm -1 ) for the 12-electron isoelectronic 
series C2, BN, and BeO obtained from a Dunham analysis of potential energy curves computed at 
the W4 and related levels. 





W4[up to 


W4 


W4 


W4 


W4 W4 6 


W4 6c 


W4.2 6c 


W4.3 C 


Expt. d 




CCSD(T) 


[up to 


[up to 


[up to 














f c,nr\ 


CCSD(T)] CCSDT] CCSDTQ] 












c 2 (X *E+) 


r e 


124.40 


124.04 


124.00 


124.19 


124.21 


124.22 


124.22 


124.22 


124.2440 


U) e 


1861.29 


1872.11 


1874.62 


1859.09 


1857.63 


1857.43 


1856.55 


1855.84 


1855.0142 


Ul e X e 


12.51 


12.48 


12.63 


13.44 


13.45 


13.45 


13.50 


13.53 


13.5547 




-0.020 


-0.019 


-0.039 


-0.081 


-0.082 


-0.082 


-0.083 


-0.088 


-0.1321 


a e 


0.01724 


0.01726 


0.01728 


0.01783 


0.01783 


0.01783 


0.01788 


0.01789 


0.018013 




BN (a X S+) 


r e 


126.85 


126.43 


127.17 


127.33 


127.38 127.42 


127.43 


127.43 


127.33 


127.45081 


U) e 


1742.85 


1755.27 


1719.62 


1713.06 


1711.08 1707.50 


1707.28 


1707.64 


1709.81 


1705.4032 


U) e X e 


16.33 


16.08 


12.29 


11.99 


11.76 11.40 


11.40 


11.37 


11.68 


10.55338 




-0.505 


-0.511 


-0.004 


-0.013 


-0.015 -0.015 


-0.015 


-0.024 


0.036 




a e 


0.01669 


0.01663 


0.01665 


0.01638 


0.01627 0.01618 


0.01617 


0.01611 


0.01632 


0.013857 




BeO (X *£+) 


r e 


133.63 


132.89 


132.85 


132.98 


132.88 133.01 


133.07 


133.07 




133.09 


Ul e 


1479.86 


1497.59 


1507.02 


1494.41 


1501.65 1493.43 


1488.55 


1488.55 




1487.32 


Ul e X e 


12.05 


12.14 


10.95 


12.12 


12.34 11.78 


11.69 


11.66 




11.83 


Oe 


0.01889 


0.01900 


0.01810 


0.01893 


0.01889 0.01881 


0.01884 


0.01882 




0.0190 



"For C2, BN, and BeO DBOC contributions are found to be insignificant and are not included in the theoretical values. 
'Diffuse functions are added to electronegative atoms (N and O) in all the valence post-CCSD(T) steps, see text. c The 
inner-shell correlation contribution is extrapolated from RCCSD(T)/aug-cc-pwCVQZ and RCCSD(T)/aug-cc-pwCV5Z 
calculations, see text. d C 2 from RefEO, BN from RefP±] an <i BeO from Re#H 
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TABLE VII: Theoretical" and experimental equilibrium geometries, harmonic frequencies, and 
fundamental frequencies (in pm, degrees, and cm" 1 ) for H2O, CO2, and O3. 





r e 8 e o)i u)2 ojz v\ V2 V3 Ref. 




CO2 


W4 
Expt. 


115.98 1351.9 670.5 2396.4 1331.5 666.4 2350.2 
116.00 1353.8 672.9 2396.5 1332.9 672.9 2349.1 E3 




H a O 


W4 

W4 6 

Expt. 

Expt. 

Expt. 


95.76 104.50 3834.6 1648.3 3945.7 3655.8 1595.3 3754.3 
3835.2 1648.0 3946.3 3656.4 1595.2 3755.0 
95.72 104.52 3832.2 1648.5 3942.5 ISO 
95.84 104.44 3831.8 1647.8 3942.4 El 

3657.1 1594.7 3755.9 EH 




o 3 


W4[up to CCSDT] 

W4[up to CCSDTQ] 

W4 

W4* c 

Expt. 

Expt. 


1153.6 717.9 1085.7 1126.5 704.9 1042.1 
1137.0 715.5 1070.2 1105.3 701.5 1020.3 
127.15 116.82 1133.9 715.2 1067.1 1100.4 701.5 1016.0 

1078.1 

127.17 116.78 1134.9 716.0 1089.2 El 
127.17 116.75 1133.3 715.0 1087.3 El 



a All the theoretical frequencies are computed at the W4-optimized geometries; unless otherwise indicated the DBOC 
contributions are not included in the theoretical values, cfr. footnote b. ^Including a DBOC correction at the 
HF/aug'-cc-pVTZ level of theory. c W4*=W4[up to CCSDT]+f 4 /cc-pVTZ(no / 1 d)+f 5 /cc-pVDZ(no d) see text. 
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TABLE VIII: Dipole moments (in a.u., 1 a.u. = 2.5416 Debye). 





W4[up to 


W4 


W4 


W4 


W4 


Expt. 


Ref. 




CCSD(T) 


[up to 


[up to 


[up to 










fc,nr] 


CCSD(T)] CCSDT] CCSDTQ] 








BH 


0.550 


0.552 


0.551 


0.551 




0.50 


I5T1 


HF 


0.709 


0.709 


0.709 


0.709 


0.709 


0.7185 a 


(ii 


HC1 


0.433 


0.430 


0.429 


0.429 


0.429 


0.4362±0.0001 


1751 


CO 


0.044 


0.046 


0.044 


0.048 


0.049 


0.0432 


(il 


cs 


0.774 


0.777 


0.771 


0.774 


0.778 


0.770 


(il 


BeO 


2.496 


2.484 


2.469 


2.470 


2.484 


2.2680 


(il 


BN 


0.796 


0.796 


0.850 


0.811 


0.799 






H 2 


0.729 


0.730 


0.730 


0.729 


0.729 


0.730±0.002 


75 


Oa 


0.218 


0.218 


0.222 


0.214 


0.213 


0.2100 


75 


MSD 6 


0.035 


0.034 


0.031 


0.031 


0.033 






MAD 6 


0.038 


0.038 


0.035 


0.035 


0.038 






RMSD 6 


0.083 


0.079 


0.074 


0.074 


0.079 






MSD C 


0.001 


0.000 


0.004 


-0.001 








MAD C 


0.004 


0.002 


0.010 


0.004 








RMSD C 


0.005 


0.002 


0.018 


0.006 









gives 0.707 a.u. Error statistics with respect to experiment. c Error statistics with respect to the W4 values. 
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TABLE IX: Static Polarizabilities (in a.u., 1 a.u. = 1.48185x 10" 25 cm 3 )". 







W4[up to 


W4 


W4 


W4 


W 4 


Expt. 


Ref. 






CCSD(T) 


[up to 


[up to 


[up to 












1 

/c,nr] 


CCSD(T)] CCSDT] CCSDTQ] 








BH 




20 87 


20.78 


20.78 


20.77 












23 14 


23.03 


23.05 


23.06 












21.63 


21.53 


21.53 


21.53 








HF 




5.15 


5.14 


5.14 


5.14 


5.14 


5.08 


b 






6.29 


6.28 


6.28 


6.29 


6.29 


6. 10 


b 






5.53 


5.52 


5.52 


5.53 


5.53 


5.52 c 


b 


HC1 




16.68 


16.65 


16.65 


16.65 


16.65 










18 32 


18.30 


18.30 


18.30 


18 30 










17.23 


17.20 


17.20 


17.20 


17.20 


17.39^ 


77 


00 




11 85 


11.81 


11.80 


11.81 


11 81 


11 86 


78 




OLzz 


15 42 


15.38 


15.36 


15.40 


15 42 


15 51 


78 






13 04 


13.00 


12.99 


13.01 


13 01 


1.3 08 


78 


OS 




23.74 


23.62 


23.63 


23.62 


23. 63 










37.93 


37.81 


37.79 


37.87 


37.85 










28.47 


28.35 


28.35 


28.37 


28.37 






A 2 




39 91 


39.66 


39.66 


39.65 


39 60 








"J" 


67.86 


67.69 


67.77 


67.52 


67.68 










49 23 


49.00 


49.03 


48.94 


19 00 






CI2 




24.86 


24.81 


24.81 


24.81 


24.81 


24.43 


79 






41.74 


41.66 


41.68 


41.65 


41.65 


44.54 


79 






30 48 


30.43 


30.44 


30.43 


30 43 


31 ll e 


79 


H2O 




9 20 


9.18 


9.18 


9.18 


9 18 










9.84 


9.83 


9.83 


9.84 


9.81 










9 48 


9.47 


9.47 


9.48 


9 18 










9 50 


9.49 


9.49 


9.50 


9 50 


9 642-^ 


77 


CO-? 




12.74 


12.74 


12.73 


12.75 


12.75 










26.55 


26.55 


26.51 


26.57 


26.61 










17 34 


17.34 


17.32 


17.36 


17 37 


17 ol 9 


77 


N2 




10.19 


10.15 


10.15 


10.16 


10.16 










14.78 


14.73 


14.76 


14.71 


14.72 










11.72 


11.68 


11.69 


11.68 


11.68 


11. 744^ 


75 


p„ 

r 2 




6.38 


6.37 


6.37 


6.37 


6.37 










12.34 


12.32 


12.36 


12.26 


12.26 










8.37 


8.35 


8.36 


8.33 


8.33 


8.38 


80 


f„ 
^■2 




36 88 


38.68 


22.86 


23.82 


23 86 








O-zz 


31.80 


31.81 


26.36 


26.33 


26.42 










35.18 


36.39 


24.03 


24.66 


24.71 






BeO 




31.46 


31.36 


30.63 


31.30 


31.59 










33.39 


33.05 


34.01 


32.91 


31.90 








a 


32.10 


31.92 


31.76 


31.84 


31.70 






o 3 


a X x 


11.75 


11.72 


11.72 


11.73 


11.73 








a yy 


31.26 


31.27 


31.03 


30.81 


30.74 








a-zz 


14.18 


14.14 


14.10 


14.11 


14.11 








a 


19.06 


19.04 


18.95 


18.88 


18.86 


21.7 


75 




MSD ; 


0.72 


0.72 


0.01 


0.01 










MAD' 


0.73 


0.75 


0.17 


0.06 










RMSD* 


2.59 


2.89 


0.47 


0.20 









a a = ^(a xx + ayy + OL zz ). 'Expt. value with zero-point vibrational corrections taken from ReflSU. c RefPH gives: q=5.4, and 
RefESgives .5=5.601 a.u. d Ref EH gives: o=17.55 a.u., Ref. Ell gives: 0:^=16.13, a zz =21.12, and a=17.75 a.u. e Ref M gives: 
3=31.1 a.u. -fRefEU gives: o=9.8 a.u. 9 Ref PI gives: (5=17.7 a.u., Ref. EH gives: o=17.88 a.u., and Ref. EH gives: o=19.64 a.u. 
h Refl EH gives: a xx = 9.79, oj zz = 16.06, and a=11.88 a.u., and RefPH gives: a=11.92 a.u. 'Error statistics with respect to the 
W4 values. 
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TABLE X: Static Hyperpolarizabilities (in a.u., 1 a.u. = 8.6392x 10~ 33 cm 4 statvolt~ 1 ) a . 







W4[up to 


W4 


W4 


W4 


W4 






CCSD(T) 


[up to 


[up to 


[up to 








fc,nr] 


CCSD(T)] CCSDT] CCSDTQ] 




BH 


ftzzz 


21.81 


22.39 


22.05 


22.16 






ftxxz 


69.95 


68.92 


68.88 


69.06 









70.85 


69.27 


69.43 


69.57 




HF 


fizzz 


8.75 


8.77 


8.79 


8.82 


8.83 




fixxz 


1.03 


1.03 


1.04 


1.05 


1.05 







6.49 


6.50 


6.52 


6.55 


6.55 


HC1 


fizzz 


9.77 


9.78 


9.81 


9.81 


9.81 




fixxz 


0.29 


0.33 


0.32 


0.31 


0.31 







6.21 


6.26 


6.28 


6.26 


6.26 


CO 


fizzz 


28.63 


28.45 


28.32 


28.67 


28.5 




ftxxz 


7.66 


7.60 


7.56 


7.64 


7.67 







26.37 


26.19 


26.06 


26.38 


26.31 


CS 


fizzz 


15.95 


15.22 


15.62 


15.90 


12.81 




fixxz 


4.79 


4.48 


4.38 


4.61 


4.72 







15.31 


14.50 


14.63 


15.07 


13.35 


H 2 I3 ZZZ 


11.86 


11.90 


11.90 


11.96 


11.96 




fixxz 


4.62 


4.62 


4.63 


4.66 


4.66 




fiyyz 


9.48 


9.51 


9.53 


9.56 


9.56 







15.58 


15.62 


15.64 


15.70 


15.70 


O3 


fizzz 


-17.32 


-17.22 


-17.40 


-17.43 


-17.44 




fixxz 


-3.61 


-3.59 


-3.61 


-3.62 


-3.63 




fiyyz 


1.66 


1.92 


-5.32 


-7.22 


-7.04 







-11.56 


-11.33 


-15.80 


-16.96 


-16.87 




MSD 6 


0.77 


0.70 


0.22 


0.18 






MAD 6 


0.77 


0.70 


0.32 


0.20 






RMSD 6 


2.26 


2.26 


0.81 


0.75 





a /3 = |(/f II2 + flyyz + Pzzz)- 6 Error statistics with respect to the W4 values. 
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-0.6 -0.3 0.0 0.3 -5 5 15 25 

(a) Ar e (pm) (b) Aw c (cm -1 ) 




-1.5 -0.5 0.5 -0.4 -0.2 0.0 0.2 0.4 

(c) Auj e x e (cm* 1 ) (d) Au e y e (cm -1 ) 



FIG. 1: Normal distribution of errors between theory and experiment for r e , u) e , uj e x e , and u> e y e . 
The centers of the Gaussians for the W4[up to CCSD(T)], W4[up to CCSDT], W4[up to CCSDTQ], 
and W4 methods, respectively, are given in parenthesis, (a) equilibrium bond lengths (-0.16, -0.17, 
-0.04, -0.03 pm), (b) harmonic frequencies (7.5, 8.1, 1.0, and 0.6 cm -1 ), (c) first-order anharmonic 
corrections (-0.30, -0.26, -0.08, -0.07 cm -1 ), and (d) second-order anharmonic corrections (-0.03, 
-0.03, -0.04, -0.04 cm- 1 ). 
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